The Lorentz singular value decomposition and its applications to pure states of 3 

qubits. 
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All mixed states of two qubits can be brought into normal form by the action of SLOCC operations 
of the kind p' = (A ® B)p{A ® . These normal forms can be obtained by considering a Lorentz 
singular value decomposition on a real parameterization of the density matrix. We show that the 
Lorentz singular values are variationally defined and give rise to entanglement monotones, with as 
a special case the concurrence. Next a necessary and sufficient criterion is conjectured for a mixed 
state to be convertible into another specific one with a non-zero probability. Finally the formalism 
of the Lorentz singular value decomposition is applied to tripartite pure states of qubits. New proofs 
are given for the existence of the GHZ- and W-class of states, and a rigorous proof for the optimal 
distillation of a GHZ-state is derived. 
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Local probabilistically reversible operations cannot af- 
fect the intrinsic nature of the entanglement present in 
a system. It is therefore interesting to look for the most 
general local operations that wash out all the local infor- 
mation such that only the non-local character remains. 
For a pure entangled state of two qubits for example, it 
is well-known that it can be locally transformed into a 
Bell state, which is indeed the only pure state for which 
the local density operators do not contain any informa- 
tion. Recently, a similar result was derived in the case 
of mixed states Q. The key ingredient of the analysis 
was the existence of a Lorentz singular value decompo- 
sition. In this report some new interesting properties 
of this Lorentz singular value decomposition are derived 
and it is shown how it is related to the existence of entan- 
glement monotones. Furthermore it leads to a criterion 
for a mixed state to be convertible into another specific 
one with a non-zero probability. It also leads to a trans- 
parent derivation of all different normal forms for pure 
states of three qubits: a pure state of three qubits is 
indeed uniquely defined, up to local operations, by the 
two qubit density operator obtained by tracing out one 
particle. It will be shown how the so-called GHZ- and W- 
type states arise. We will also give a rigorous proof 
of the optimal way of distilling a GHZ-state, confirming 
the results of Acin et al. llloll . 



I. THE LORENTZ SINGULAR VALUE 
DECOMPOSITION 

Let us consider a mixed state of two qubits and inves- 
tigate the orbit generated by probabilistically reversible 



SLOCC operations of the kind 



(1) 



where A, B are complex 2x2 matrices of determinant 1 
and p' is unnormalized. It will turn out very convenient 
to work in the real i?-picturc defined as 



1 ^ 



ij=0 

with {ai} the Pauli spin matrices. As shown in the 
determinant 1 SLOCC operations (|^) in the p-picture be- 
come proper orthochronous Lorentz transformations in 
the i?-picture: 

R' = Lj\RL^ 

Indeed, it is a well known accident that SL{2,C) ~ 
50(3, 1). Note that local operations are very transparent 
in the i?-picture: operations by Alice amount to opera- 
tions on the row space of R, i.e. left matrix multipli- 
cation, while operations by Bob result in column opera- 
tions. 

Let us now state a more refined version of the central 
theorem of 

Theorem 1 The 4x4 matrix R with elements R^ = 



Tr [pui ® <7j) can be decomposed 



R — L/i YjL/n 



with Li, L2 finite proper orthochronous Lorentz transfor- 
mations, and S either of unique real diagonal form 



^In comparison with the original theorem we introduced a 
more refined classification in the case of non-diagonalizable 
R. 
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/ So 



V 



Si 



S2 



S3 



with So > si > S2 > jsal and S3 positive or negative, or 
of the form 



S = 



/ 



V 



with unique a, b, c, d obeying one of the four follow- 
ing relations: [b — c = a/2); ((d = = c) A {b = a)); 
{{d = 0^b)A{c = a)); ((d 0) A (a = = c)). 

It is interesting to note that the Lorentz singular val- 
ues are the only invariants of a state under the SLOCC 
operations (|^). 

The diagonalizable case is generic, and a diagonal R 
corresponds to a Bell-diagonal state. The existence of 
the non-diagonal normal forms is a consequence of the 
fact that the Lorentz group is not compact: these non- 
diagonal normal forms can only be brought into diagonal 
form by infinite Lorentz transformations. Nevertheless, 
even in those cases the Lorentz singular values are well 
defined and given by: 



[so, si, S2, S3] = [v/(a - 6) (a - c), ^(a - b){a - c), d, -d] 

The corresponding normal form of the non-diagonalizable 
case in the p-picture is given by: 



filtering operations by one party is enough to distill a sin- 
glet out of an non-maximally entangled pure state: Alice 
or Bob can apply the filter corresponding to the Lorentz 
transformation given by the inverse of R. 

The success of the ordinary singular value decompo- 
sition is to a large extent the consequence of the nice 
variational properties of the singular values: the sum of 
the n largest singular values is equal to the maximal in- 
ner product of the matrix with whatever n orthonormal 
vectors. Interestingly, a similar property holds for the 
Lorentz singular values: 

Theorem 2 The Lorentz singular values sq > si > S2 > 
I S3 1 of a density operator R are variationally defined as: 



/I 



So = mill Tr L1RL2 

L\ ,L2 



So — Si = min Tr 



LIRL2 



So — Si — S2 = min Tr 

Li,L2 



( /I 

V 



So — si — S2 -I- S3 = min Tr 

Li ,L2 

where Li,L2 are proper orthochronous Lorentz transfor- 
mations. 



a — b d 
d a ^ c 



Proof: We will give a proof for the fourth identity and 
the other proofs follow in a completely analogous way. 
An arbitrary Lorentz transformation can be written as 



The four distinct non-diagonal normal forms correspond 
to the following states: 

• {b = c = a/2): these are rank 3 states (rank 2 iff 
{d = b = c)) with the strange property that their 
entanglement cannot be increased by any global 
unitary operation 

• ((d = = c) A (6 = a)); ((d = = &) A (c = a)): p 
is separable and a tensor product of the projector 
diag[l; 0] and the identity. 

• ((d = 0) A (a = 6 = c)): p is the separable pure 
state diag[l; 0; 0; 0]. 

Let us now consider the case of a generic pure state for 
which we always have the relation sq = si = S2 = —S3. 
This implies that R itself is a Lorentz transformation, 
up to a constant factor that will turn out to be the con- 
currence; the singlet state for example is represented in 
the i?-picture by i? = diag[l; 1; 1; —1]. This clarifies why 



L = 



V 



/ cosh(a) sinh(a) 
sinh(Q!) cosh(a) 



w 



1 



where V and W are orthogonal 3x3 matrices of determi- 
nant 1. There is no restriction in letting R be in normal 
diagonal form, and therefore we have to find the mini- 
mum of 



Tr 



cosh(Q!) sinh(Q!) 
sinh(a) cosh(a) 



1 



W 



V 



over all V, W, a. Using the variational properties of 
the ordinary singular value decomposition and the fact 
that the Lorentz singular values are ordered, it is im- 
mediately clear that an optimal solution will consist in 
choosing W = I3, V = diag[— 1; — 1; 1] and a = as 
cosh(a) > sinh(a) and sq > si. This ends the proof. □ 
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II. LOCAL INVARIANTS VERSUS 
ENTANGLEMENT 



Let us now investigate how the local invariant Lorentz 
singular values are related to the concept of entangle- 
ment. Inspired by theorem we define the quanti- 
ties Mi{p) = max(0,— (so — si — S2)) and A'hip) = 
max(0, — (so — si — S2 -\- S3)). As they are solely a func- 
tion of the non-local invariants of the density operator, 
we suspect them to be related to the amount of entan- 
glement present in the considered state: 



Theorem 3 Ml (p) = max(0, - 
M2{p) = max(0, -(sq - si - S2 ^ 
monotones. 



'{sq — Si — S2)) and 
S3)) are entanglement 



Proof: A quantity M{p) is an entanglement monotone 
iff its expected value decreases under the action of 
every local operation. Due to the variational characteri- 
zation of the quantities (sq — Si — S2) and (sq — si — S2-I-S3), 
it is immediately clear that both Mi and M2 are decreas- 
ing under the action of mixing. It is therefore sufficient 
to show that for every A < I2, A ~ \/ 12 — A, it holds 
that 



M,{p) > 

Tr {{A (g) I)p{A (g> ly) M, 
+Tr {{A (g) I)p{A (g> /)t) M, 



{A(g)I)p{A(g,I)^ 
TT{{A(g>I)p{A(g)iy)^ 

( {A®i)p{A®iy 

i Tr((i(g)/)p(i(g)/)t 



It should be clear from the previous discussion that the 
following identity holds: 



det(A)M,(p) 



(A /)p(A (g) /)t 
Tr {{A ® I)p{A ® /)t) ; Tr {{A ® I)p{A ® /)t) 



Indeed, A/ yJdei{A) corresponds to a Lorentz transfor- 
mations which cannot change the Lorentz singular values. 
We therefore only have to prove that 1 > det(A)-|-det(A). 
Given the singular values cri,(T2 of A, this inequality is 
1 > tTi(T2 -I- (1 — (Ti)(l — (T2) which is trivially fulfilled. □ 
Both Ml and M2 , linear functions of the Lorentz sin- 
gular values, are therefore entanglement monotones that 
are analytically calculable for mixed states: we did not 
use the concept of convex roof formalism. It turns out 
that M2 is equivalent to the concurrence of a state as 
introduced by Wootters || : 



C{p) 



0,-2(^0 - |si| 



IS2I + S3) ) = -M2{p) 



There is indeed a strong relation between the Lorentz 
singular values and the eigenvalues {A^} of the operator 
\/ {(7y (g (Ty)p'^{(jy ® ay)p Introduccd by Wootters |^ : 







/Ai\ 




A2 


) 


A3 







Together with the negativity the above entangle- 

ment monotones are the only ones for which an analytical 
expression exists for whatever mixed two-qubit state. 

The existence of entanglement monotones is interest- 
ing as it gives necessary conditions for one state to be 
convertible into another one by LOCC operations with 
probability 1. It is still an open problem to find the suffi- 
cient conditions for the convertibility of one mixed state 
into another one, although this was solved for pure states 
[^"ll . If we relax the constraints that the conversion has 
to succeed with unit probability, the above formalism 
can give us some answers in the case of mixed states. 
We have indeed shown that a generic state can always 
be brought into Bell-diagonal form by the SLOCC oper- 
ations (|l|). The problem of one state to be convertible 
into another one with a non-zero probability is therefore 
reduced to the question whether one Bell diagonal state 
can be transformed into another one. A Bell-diagonal 
state is uniquely defined under the SLOCC operations 
(|l|), and therefore the only local tool remaining is mix- 
ing. Numerical and theoretical investigations indicate 
that a given Bell diagonal state can only be converted 
into another one iff this last one is a mixture of the origi- 
nal Bell-diagonal state with a separable state, although a 
general proof has not been found. We conjecture however 
that this is always true: 

Conjecture 1 A two-qubit state pi can probabilistically 
be converted into the state p2 iff the Bell-diagonal normal 
form of p2 is a convex sum of a separable state and the 
Bell-diagonal normal form of pi. 

It is clear that a trivial procedure exists to implement 
this conversion with unit efficiency: mix the state with 
one that can be locally made. Let us for example inves- 
tigate whether the Bell-diagonal pi with ordered eigen- 
values {Ai} can be transformed into the Bell-diagonal p2 
with ordered eigenvalues {pi}. We can restrict ourselves 
to mixing with separable Bell diagonal states lying on 
the boundary of the entangled and separable states, and 
these have their largest eigenvalue equal to 1/2. Un- 
der the assumption of our conjecture, conversion is pos- 
sible iff the following constrained system of equations in 
X, y, z, t, P has a solution: 




= (1 




/1/2 

y 

z 

t 



(2/,z,i>0) {v + z + t = l/2) 



where P3 is a 3 x 3 permutation matrix. This system can 
readily be solved. Not surprisingly, there is a close rela- 
tion between majorization and the above set of equations. 
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Note also that a pure entangled state can be converted 
probabilistically into whatever mixed state. 



III. PURE STATES OF THREE QUBITS 

Let us now apply the Lorentz singular value decom- 
position on the problem of determining the different 
classes under SLOCC operations of pure three qubit 
states. We will not consider the states that have a ten- 
sor product structure, as these are not truly tripartite. 
Therefore we know that taking the partial trace over 
whatever party will result in a rank 2 density opera- 
tor of two qubits. Due to theorem ^, we know that 
this density operator can be brought into one of two 
normal forms by SLOCC operations of two parties: a 
Bell diagonal pi ~ p\tjA){ip^\ -|- (1 — with 
1^-+) = (|00) + |ll))/v/(2);|V'-) - (|00) - |ll))/v/(2); 
or a quasi distillable p2 = \(t)+){(j)+\ + |00)(00| with 
|(/,+ ) = (|01) + |10))/vT2), but this last case is clearly 
not generic. 

Purification of this second normal form directly leads 
to the normal form |100) + |010) + |001), called the W- 
state As shown in j^], this state has maximal two- 
partite entanglement distributed over all parties. 

Purification of the Bell diagonal state results in 

IV') = a(|000) + |110)) + 6(|001) - 



If the third party now applies the local operation A = 
a b 



^ _^ ] , the GHZ state |000) + |111> with maximal 

true tripartite entanglement is obtained. 

Therefore we have given an alternative proof of the 
following theorem by Diir, Vidal and Cirac 

Theorem 4 Every pure tripartite entangled state can 
he transformed to either the GHZ, or the W-state, by 
SLOCC operations. 

Note that the SLOCC operations bringing a generic 
pure state to the GHZ form are not unique but con- 
sist of a four-parameter family. This happens because 
a pure tripartite state has 14 degrees of freedom and the 
three Lorentz transformations have 18 independent pa- 
rameters. Indeed, \i A® B ® CIV') = \GHZ), then also 



a rigorous proof of the same result is presented in the 
appendix . 

A similar non-uniqueness exists in the case of distilling 
a state of the class of W-states to the W-state. Indeed, 
\i A® B ® C\'4>) — \ W), then the most general symmetry 
operations are given by 

A' ®B' ®C'\i}) = \W) 



A' 
B' 



X y 
l/x 

X z 



A 
B 



l/x 

^1 _ I X -x{y + z) 
^ - ' l/x 



C 



with x,y, z arbitrary complex numbers. As every ma- 
trix can be written as the product of a unitary matrix 
and an upper triangular matrix (this is the so-called QR- 
decomposition) , there are enough degrees of freedom to 
make whatever one out of A' ,B' or C" equal to a unitary 
matrix. Numerical investigations reveal that one of these 
three possibilities is also the optimal choice in the sense 
that it will yield a distillation protocol that produces the 
W-state with the highest possible probability. Therefore 
the optimal distillation protocol of a W-state consists of 
two parties applying a local filtering operation, while one 
party performs a local unitary operation. 

Finally, a natural question arises as how the previous 
results generalize to the case of mixed states. Due to the 
fact that the rank of a density matrix corresponding to 
a mixed state is higher then 1, it is immediately clear 
that no SLOCC operations can exist that yield a rank 1 
GHZ-state. In [^| it is shown that the optimal SLOCC 
operations in the case of mixed states are those that pro- 
duce a unique state from a given state such that all its 
local density matrices are equal to the identity. Note that 
the GHZ state is the only pure state with this property 
in the 3-qubit case. 

A second question concerns the generalization of the 
class of pure W-states to mixed states: is the W-class of 
mixed states of measure zero? This question was solved 
in 1 13 1 (see also |14 for a simple derivation and general- 
izations), where it was shown that the W-class of mixed 
state is not of measure zero. 





1/a 



A( 





1/6 



B< 



1/ab 






ab 



C\iIj) = \GHZ) 
(2) 



with a and b complex numbers. The single-copy distilla- 
tion of a GHZ state is therefore not unique. The prob- 
ability by which an SLOCC operation (|l) produces the 
desired result can therefore be optimized such as to yield 
the optimal single-copy distillation protocol. This opti- 
mal procedure was previously found by Acin et al. pO[; 



APPENDIX A: OPTIMAL DISTILLATION OF 
THE GHZ STATE 

The most general local procedure of distilling a GHZ- 
state out of a single copy of a pure state consists of a 
multi-branch protocol in which different branches con- 
sist of different SLOCC operations connected through 
equation (|). There is no restriction in taking all 
{Ai}, {Bi}, {Ci} to have determinant 1, and the SLOCC 
operations corresponding to each branch are of the form 
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q,T^'^\GHZ) 

ai 
l/a, 

h 
1/6, 

l/aih 

Qih 



Here r is the 3-tangle [|11|,|_2| of V' and qi is a real pro- 
portionality factor such as to assure that all the branches 
together are implementable as a part of a POVM. This 
leads to a necessary (but generally not sufficient) condi- 
tion: 



B B,®C a,<h 



(Al) 



Each branch yields the GHZ-state with probability 




and therefore the total probability is given by ^Jt ^ ^- q^ 
which has to be maximized. Due to the condition ( [Al| ), 
an upper bound on this probability can readily be de- 
rived. It will turn out that this upper bound is achiev- 
able by a 1-branch protocol. Defining pi ~ Qi/{J2iQi)j 
it holds that the total probability is bounded by 



Therefore only the (real and positive) diagonal ele- 
ments of ZqUPuWzI are of importance and let us 
write them in the vector zq. Similarly, we write 
the diagonal elements of Di in the vector di = 
[l;|aAP;l/|6,p;|a,|2;l/|a,|2;|5,|2;l/|aAP;l], and the 
extremal relations become: 



: = df diag[0, 1, 0, 1, -1, 0, -1, 0]zo 
= diag[0, 1, -1, 0, 0, 1, -1, 0]zo 
dfzo 



(A3) 



where /i is the Lagrange multiplier corresponding to the 
condition J2iiVPi^^ ~ ^- "^^^^ forms sets of each time 3 
equations for 2 unknowns 0^,6^, which can be shown to 
have exactly one solution. Indeed, the first and second 
equation lead to 



zo{5) + z„{7)/\b. 



^o(4) 
^o(3) - 



-zo{2)\b,\^ 
^o(7)/|a,p 



zo(6) + zo(2)|a, 



(A4) 



max 



(A2) 



where Xyna^iX) denotes the largest eigenvalue of opera- 
tor X. An upper bound is therefore obtained by mini- 
mizing this largest eigenvalue. Therefore the standard 
techniques for differentiating the eigenvalues of a ma- 
trix have to be used [ p^ : given a Hermitian matrix X, 
its eigenvalue decomposition X = UEW and its varia- 
tion X, then the variation on its eigenvalues is given by 
E = di&g{U^XU}. Here we take 



X 



Zo 
D, 



Ao®Bo® Co 



c|2 



Note that varying the free parameters {ai,hi,pi} only 
affects D and not Zq. In the case of an extremal max- 
imal eigenvalue all variations Amax = Tr {ePi-^ with 
Pii — diag[l; 0; 0; 0; 0; 0; 0; 0] have to be equal to zero: 

Tr(((5i?)ZoC/Piit/t4) =0 

The following identities are easily verified: 
3D 2 

— = -diag[0, 1,0, 1,-1, 0,-1,0] A 
SD 2 

— = -diag[0, 1,-1, 0,0, 1,-1,0] A 
obi b. 



SD 

SJWi 



Let us analyze how these equations behave. When 
6i — > th en t he solution of the first equation goes like 
ki| ^ '^l\fWi\ and when a; then ~ l/|aip. Ex- 
actly the opposite happens in the case of the second equa- 
tion, and due to this different asymptotic behaviour it is 
assured that both curves cross and therefore at least one 
solution exists for all (real positive) values of zq. More- 
over there is always at most one solution. To prove this, 
we first note that \ai\ and can be scaled such that 
both curves cross at the value (1,1), and we call these 
rescaled variables (x,y) and zp . The hyperbola xy = \ 
crosses both rescaled curves ( |A4D at (1, 1). Moreover it 
is trivial to check that the hyperbola does not cross any 
of the rescaled curves anymore in the first quadrant (this 
amounts to solving a quadratic equation) , and due to the 
asymptotic behaviour one curve lies below and the other 
one above the hyperbola (except in (1,1)). Therefore 
both rescaled curves have exactly one crossing. There- 
fore for all (real positive) values in zq, there is always 
exactly one real solution for |ai|, and as is inde- 
pendent of the index ?, all \ai \ are equal to each other and 
the same applies to the Therefore at most the phase 
of the constants {ai,6i} varies in different branches, and 
as this amounts to local unitary operations we conclude 
that all branches are equivalent and can be implemented 
by a one-branch protocol. This implies that the upper 
bound (A2) can be reached. 

In the case of a one branch protocol, the eigenvectors of 
X can be calculated analytically as X becomes a tensor 
product of 2 X 2 matrices. Given particular determinant 
1 transformations A, B, C and taking a, b to be real, the 
eigenvector v corresponding to the largest eigenvalue of 



5 



a 
1/a 



the matrix YY'' with Y - 

^ir^ abl^ happens to be v 



b 

l/b 

1 vs with 



Acin, E. Jane and G. Vidal. FV is very grateful to Hideo 
Mabuchi and the Institute of Quantum Information at 
Cahech, where part of this work has been done. 



ai 
a2 



2y/AiiA22 - 1 



2^ B11B22 
as = 2\/ C11C22 



1 



(3i = Ana^ - ^22/0^ 
/32 = Bnfo' - B22/b^ 
Pz^Cii/iabf -C22{abf 



The conditions (A3) then imply that ft/ai is a constant 
for all i — 1..3: 

Aiia2-A22/a2 Bnb^-B22/b^ Cn / (ab)^ - C 22 (ab)^ 



VA11A22 - 1 



V B11B22 — 1 



vCnC; 



22 



1 



These two equations have to be solved in the unknowns 
a and b. b can readily be written in function of a through 
one of those, and then a sixth order equation in the re- 
maining unknown results. As shown above, only one 
solution corresponding to a physical solution for a and 
b exists, and this solution can easily be solved numer- 
ically. The optimal local filtering operations and the 
maximal probability of making a GHZ-state (an entan- 
glement monotone Q) can then easily be calculated. 
The solution obtained is completely equivalent to the one 
of Acin et al. jl^, although their proof did not include 
the uniqueness of the solution and needed exhaustive nu- 
merical calculations. 

Note that the procedure outlined here is equally applica- 
ble to the problem of distilling a GHZ-state in a higher 
dimensional system. 
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